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ABSTRACT 

We present SubMillimeter- Array observations of a Keplerian disk around the Class 
I protobinary system L1551 NE in 335 GHz continuum emission and submillimeter line 
emission in 13 C0 (J =3-2) and C 18 (J =3-2) at a resolution of ~120 x 80 AU. The 
335-GHz dust-continuum image shows a strong central peak closely coincident with 
the binary protostars and likely corresponding to circumstellar disks, surrounded by a 
~600 x 300 AU feature elongated approximately perpendicular to the [Fe II] jet from 
the southern protostellar component suggestive of a circumbinary disk. The 13 CO and 
C 18 images confirm that the circumbinary continuum feature is indeed a rotating 
disk; furthermore, the C 18 channel maps can be well modeled by a geometrically-thin 
disk exhibiting Keplerian rotation. We estimate a mass for the circumbinary disk of 
~0. 03-0. 12 Mq, compared with an enclosed mass of ~0.8 M & that is dominated by the 
protobinary system. Compared with several other Class I protostars known to exhibit 
Keplerian disks, L1551 NE has the lowest bolometric temperature (~91 K), highest 
envelope mass (~0.39 Mq), and the lowest ratio in stellar mass to envelope + disk + 
stellar mass (~0.65). L1551 NE may therefore be the youngest protostellar object so far 
found to exhibit a Keplerian disk. Our observations present firm evidence that Keplerian 
disks around binary protostellar systems, "Keplerian circumbinary disks", can exist. We 
speculate that tidal effects from binary companions could transport angular momenta 
toward the inner edge of the circumbinary disk and create the Keplerian circumbinary 
disk. 
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1. Introduction 

Previous interferometric observations of millimeter and submillimeter continuum emission 



in extended f>2000 AU) enveloDes (Loonev et al. 200C; 


Harvev et al. 2003a: 


Harvev et al. 2003b: 


Enoch et al. 2009: 


Enoch et al. 2011 


: Joreensen et al. 2009 


■ 


Maurv et al. 20101). Sub-arcsecond- 



resolution continuum observations that resolve the compact components re veal flattened mor- 



phologies common l y attributed to circumstellar or possib ly circumbinary disks (jLoonev et al. 20001 : 



Brinch et al. 20071 : lEnoch et al. 20091 : lEnoch et al. 20 111 ) . The masses of the putative circumstellar 



disks around Class and I protostars are on the order of 0.01 - 1.0 Mq, mu ch higher than the typical 



masses of circumstellar disk s around Class II (T-Tauri) stars (~0.005 Mq) ([Andrews Williams 20051 : 



Andrews &: Williams 20071 ) . The sizes of the circumstellar di s ks (~600 AU) around all these classes 



of stars, however, are comparable (jGuilloteau Dutrev 19981 : iGuilloteau et al. 19991 : ISchaefer et al. 20 09) 



Previous interferometric observations in molecular lines have shown an i nfalling motion and /or 
little rotating mo t ion in circumstellar disks aroun d Class protostars (jTakakuwa et al. 20071 : 
Brinch et al. 20091 : lYen et al. 201G : lYen et al. 20 111 ) . Such " pseudo-disks" have been predicted by 
theoretical simulations to form in the inner regions of envelopes around protosta rs ( Mellon &; Li 20081; 
Mellon Li 20091 : IVorobvov 20091 : iMachida et al. 201ol : iMachida et al. 2011al : iMachida et al. 201 lbh 
On the other hand, interferometric observations of T-Tauri stars in millimeter molecular lines have 



revealed Keplerian disks around Class II stars (jGuilloteau &: Dutrev 199a : IGuilloteau et al. 19991 : 
Pietu et al. 20031 : ISchaefer et al. 20091 ) . These studies suggest that the transition from pseudo-disks 
to Keplerian disks occurs sometime between the Class and II stages. 

Interferometric observations of molecular lines with, especially, the SubMillimeter Array (SMA), 
which can (better) image higher transitions that trace gas at higher densities and/or temperatures, 
have been increasingly successful at finding Keplerian disks around Class I protostars. Brinch et al. 
(2007) have identified a (radius) r ~200 AU-scale Keplerian disk around L1489 IRS with a central 
stellar mass of 1.35 Mq. Observations by Lommen et al. (2008) have also revealed a r ~100 AU- 
scale Keplerian disk with a central stellar mass of 0.37 Mq around IRS 63, and a r ~200 AU-scale 
Keplerian disk with a central stellar mass of 2.5 Mq around Elias 29. Observations by J0rgensen 
et al. (2009) have shown a r ~140 AU-scale Keplerian disk around the Class I protostar IRS 43 
with a central stellar mass of 1.0 Mq. The masses of these Keplerian disks range from ~0.004 to 
~0.06 Mq, slightly larger than that around T-Tauri stars, whereas their radii are comparable to 
those of the Keplerian disks observed around T Tauri stars. 

Current theoretical simulations have trouble transforming pseudo-disks to Keplerian disks 
around protostars. Mellon & Li (2008; 2009) have performed ideal magnetohydrodynamic (MHD) 
simulations of self-similar, collapsing envelopes to investigate formation of Keplerian disks. They 
found that, for typical magnetic field strengths in interstellar clouds, magnetic braking is so efficient 
that centrifugally-supported disks cannot form. Machida et al. (2010; 2011a; 2011b) have conducted 
3-D MHD simulations of a collapsing Bonnor-Ebert sphere, including non-ideal MHD effects and 
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finite spherical radii (R c = 6.2 x 10 3 AU for a mass of 1.0 Mq). Their models show that the 
Ohmic dissipation and the depletion of the surrounding infalling envelope that anchors the magnetic 
field significantly reduce the magnetic braking and proceed the disk formation. The predicted 
masses of the disks, however, are ~a few to ~a few hundred times higher than the masses of the 
central protostar, and hence the disks are subject to further fragmentation. Thus further long-term 
simulations are required to reproduce the observed Keplerian disks around the Class I protostars. 
To provide observational insights to theoretical models, it is important to identify a Keplerian disk 
just after the transition from the pseudo-disk, and to compare its properties to those of previously- 
identified Keplerian disks and the theoretical models. 

Millimeter continuum observations have also revealed evidence for disks around binary proto 



stars, "circumbinary di s ks" (ILoonev et al 



Lim fc Takakuwa 20061 : iMaurv et al. 2010l : 



997l:lLoonev et al. 2000l : lMoriartv-Schieven et al. 2000l : 
Chen et al. 20121 ). It is intriguing to investigate whether 



Keplerian rotation is present in those circumbinary disks as well as disks around single protostars. 
Although the Class I protostars around which Keplerian disks have been found are thought to be 
single, it is possible that many and perhaps all are binary (or multiple) given the high frequency of 

binarity around pre- MS and MS solar-mass stars (jDuquennov fe Mayor 199ll; iDelfosse et al. 19991; 

Mathieu et al. 200oh . as has been suggested for L1489 IRS (jHogerheiide fe Sandell 200ol : brinch et al. 20071 1 
It is therefore important to compare the properties of circumbinary disks around known closely- 
separated binary protostellar systems with the presumed circumstellar disks of apparently single 
protostars. 

In this paper, we report observations of the Class I protostar L1551 NE in the 335 GHz 
continuum emission, along with the 13 CO ( J=3-2) and C 18 ( J=3-2) lines, at sub-arcsecond res- 
olution. Compared to the abovementioned Class I sources associated with the Keplerian disks, 
LI 5 51 NE is likely the youngest Class I object, since the bolometric temperature (T^ = 91 
K) (jFroebrich 20051 ) is more than a factor 3 lower and the amount of the surrounding envelope 
around L1551 NE is ^ one order of magnitude higher than that around the oth er protostars 



(jYokogawa et al. 20031 : iMoriarty-Schieven et al. 20061 : iTakakuwa fe Kamazaki 20 111 ). L1551 NE 



contains two 3.6-cm radio continuum sources with a projected separation of ~70 AU at a posi- 
tion angle of ~300°: the south-eastern source is referred to as "Source A" and the north-western 



source "Source B" (jReipurth et al. 2 002). Near-Infrared observations of L1551 NE have revealed 
that Source A drive s collimated [Fe II] je t s along the north-east to south-west direct ion at a posi- 
tion angle of ~60° (jReipurth et al. 200ol : iReipurth et al. 20021 : lHavashi fe Pvo 20091). and Source 
B is located at the origin of an extended 2000 AU) NIR reflection nebula (jReipurth et al. 200ol : 
Reipurth et al. 20021 : lHavashi fe Pvo 20091 ). There is little doubt that LI 551 NE is a binary pro- 
tostellar system. Our observations of this nearby (d = 140 pc; Elias 1978) object reveal for the 
first time a Keplerian circumbinary disk. In the following, we shall describe observations (§2), con- 
tinuum and molecular-line results (§3), and the analysis to identify a Keplerian disk (§4). In the 
last section (§5), we will compare the properties of the Keplerian disk around L1551 NE to those 
around the other protostars and discuss the evolution and the formation mechanism of circumstellar 
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/ Keplerian circumbinary disks. 



2. SMA Observations 

We observed L1551 NE with the extended configuration of the SubMillimeter Array (SMA) 
on 2010 January 15 and 18. The 13 CO (J=3-2; 330.587965 GHz) and C 18 (J=3-2; 329.3305453 
GHz) lines, and the continuum emission at both the upper and lower sidebands, were measured 
simultaneously. The spectral windows ("chunks") of the SMA correlator assigned to the 13 CO and 
C 18 lines had a resolution of 406.25 kHz, corresponding to a velocity resolution of ~0.37 km s _1 . 
In each sideband there are totally 48 chunks, providing a total bandwidth of 4 GHz. All the chunks 
at both sidebands except for those assigned for the 13 CO and C 18 lines were combined to make 
a single continuum channel. The frequency coverages of the continuum observation were 327.16 
- 329.24 GHz (A = 0.911 - 0.916 mm), 329.34 - 330.47 GHz (0.907 - 0.910 mm), 330.57 - 331.13 
GHz (0.905 - 0.907 mm), and 339.16 - 343.13 GHz (0.874 - 0.884 mm). The total bandwidth of the 
continuum observation was 7.74 GHz. In the following, we refer to the continuum emission taken 
with the present SMA observations as "335 GHz continuum emission" , since the mean frequency of 
the combined continuum data is ~335.4 GHz. Details of the SMA are described by Ho et al. (2004). 
Table 1 summarizes the observational parameters. The minimum projected baseline length was ~29 
kX at the C 18 frequency, and for a Gaussian emission distribution w ith a FWHM of ^7" ( ~1000 



AU), the peak flux recovered is ~10% of the peak flux of the Gaussian (jWilner &; Welch 19941 ). The 



estimated uncertainty in the absolute flux calibration is ~30%. The raw visi bility data were cal l 



brated and flagged with MIR, which is an IDL-based data reduction package (jScoville et al. 19931 ) 



The calibra ted visibility data were Fourier-transformed and CLEANed with MIRIAD to produce 



the images (jSault et al. 19951 ) 



3. Results 

3.1. 335 GHz Continuum Emission 

Figure [U shows the 335 GHz continuum image of L1551 NE at an angular resolution of / .'80 
x / /54 (P.A. = -87°), corresponding to a spatial resolution of 110 x 80 AU. The submillimeter 
continuum emission exhibits three distinct peaks. The bright central peak is located ~30 AU SE 
from Source A, and exhibits a marginal extension towards Source B (note that the central peak is 
elongated in a direction that is different to that of the synthesized beam). This central component 
appears to be closely coincident with the protobinary. From the 2-dimensional Gaussian fitting, the 
total flux density of this component is measured to be ~0.73 Jy and the deconvolved size (//90 x 



1 The SMA is a joint project between the Smithsonian Astrophysical Observatory and the Academia Sinica Institute 
of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. 
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C//49 (P. A. = -72°). The other two peaks, one to the north and the other to the south, are located 
almost symmetrically with respect to the central component. These peaks exhibit tilted [/-shaped 
morphologies (see dashed curves in Figure [1]) , which trace a (diameter) D ~600 x 300 AU-scale 
inclined ring- like feature. The position angle of this feature (~160°) is approximately perpendicular 
to the axis of the [Fe II] jets driven by Source A (P. A. = 64°; see blue and red arrows in Figure [1]). 
The total flux density integrated over the entire emission region including these three components 
is ~1.65 Jy. A similar three-peak feature, albeit at a lower S/N ratio, is also seen in the 1.3-mm 



continuum image of LI 551 NE taken with the OVRO (lMoriartv-Schieven et al. 20001 ). Without the 
benefit of observations in molecular lines, however, the interpretation made by Moriarty-Schieven 
et al. (2000) for the continuum structure is very different from ours as described in Section 4. 

Single-dish observations of L1551 NE in the 8 50 um continuum show an ext ended (= 9'/8 x 



8" 8) envelope with a total flux density of 2.78 Jy (jMoriarty-Schieven et al. 20061 ). To assess the 
contribution from the extended envelope to the more compact emission observed with the SMA, 
we created a 9 / .'8 x 8'/8 Gaussian image with a total flux density of 2.78 Jy, and sampled the 
emission with the uv-coverage of our SMA observations. We found that the envelope is almost 
entirely resolved out: we recovered a flux density of only ~7.6 mJy, less than 0.5% of the total 
flux density of the more compact emission observed with the SMA. Thus, the contribution from 
the extended envelope to the more compact structure observed with the SMA is likely negligible. 
Similarly, J0rgensen et al. (2009) found that the contribution from extended envelope emission 
measured at mm/sub- mm wavelengths in single-dish observations to that in SMA observations 
having baselines >50 kX is at most 8% for various envelope models. We estimated a total (gas + 
dust) mass (=M g+ d) for the feature detected in the 335-GHz continuum as shown in Figure 1 of 

where v is the frequency, S u the flux density, d the distance, B v {Td) the Planck function, and 
Td the dust temperature, and k v the dust opacity per unit gas + dust mass on the assumption 
that the gas-to-dust mass ratio is 100. The frequency dependence of k u can be expressed as n v = 
Kv (v / , where (3 denotes the dust-opacity index. By modeling the spectral energy distribution 
(SED) of the envelope as measured in single-dish observations at wavelengths from 12 /im to 2 
mm, Moriarty-Schieven et al. (1994) derived (3 ~1 and Barsony & Chandle r (1993) and Moriarty - 



Schieven et al. (1994) T d K. If we adopt k vo = k 12 ooghz = 0.1 cm 2 g" 1 (|Beckwith et al. 1990l ). 
then from Eq. (1) we find M g+( i ~0.047 Mq. This value is comparable with that estimated by 
Moriarty-Schieven et al. (2000) from their OVRO 1.3-mm image (~0.060 Mq) based on the same 
parameters but after correcting for their different adopted distance of ~160 pc. J0rgensen et al. 
(2009) performed a survey of low-mass protostars with the SMA and adopted a dust opacity of 
of Ossenkopf & Henning (1994) for grains with thin ice mantles coagulated at a density of n# 2 
~10 6 cm" 3 ; i.e., K850/xm = 0.0182 cm 2 g _1 and /3 = 1.7. To permit a more direct comparison with 
the results of J0rgensen et al. (2009), we adopt this dust opacity law and = 42 K, and from 
eq.(l) we find M g+C i ~0.078 Mq. J0rgensen et al. (2009) referred to simulations by Visser et al. 
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(2009) that showed dust temperatures at r=200 AU ranging from ~20 to 50 K. This range in dust 
temperatures gives M g+d = 0.063 M & (T d = 50 K) - 0.21 Mq (T d = 20 K) with the same dust 
opacity law of Ossenkopf & Henning (1994). These estimates show that M g+d is probably within 
the range of ~0.05 - 0.2 M . The same methods and the flux density of the central component 
derived from the 2-dimensional Gaussian fitting (~0.73 Jy) give the mass of the central component 
of ~0.02 - 0.09 Mq. Subtraction of this mass from M g+d gives the mass of the rest of the emission 
components of ~0.03 - 0.12 Mq. 



3.2. 13 CO (3 2) and C ls O (3 2) 

Figure [2] shows integrated intensity maps in 13 CO (3-2) and C 18 (3-2). In both 13 CO and 
C 18 0, the emission is elongated along the north-west to south-east direction. This elongation is 
approximately perpendicular to the axis of the [Fe II] jets driven by Source A, just like the ring-like 
structure seen in the continuum. Furthermore, just like this structure, the emission in both lines 
has outer dimensions of ~600 x 300 AU. Although the overall sizes and position angles of the 
major axes of molecular and ring-like continuum structures are similar, their internal structures 
are quite different. The three peaks seen in the continuum emission are not evident in 13 CO and 
C 18 0. Instead, the 13 CO emission exhibits a single peak between Source A and Source B (~40 AU 
NW from Source A), whereas the bright central peak in the continuum emission is located ~30 
AU SE from Source A. The positional accuracy of the 13 CO emission peak can be approximately 
estimated as ~beam size / signal-to-noise ratio <5 AU, and that of the continuum emission peak 
<1 AU. Thus, the positional offsets of the 13 CO emission peak from the continuum peak and the 
position of Source A are statistically significant. By contrast, the C 18 emission does not show a 
central peak, but instead shows two peaks at ~80 AU NW and ~170 AU SE from Source A, closer 
in than the two outer continuum peaks which are located ~210 AU NW and SE from Source A. The 
positional accuracy of these two C 18 emission peaks is estimated to be <13 AU, and that of the 
two continuum peaks <5 AU. Therefore, the offsets of the two C 18 emission peaks from the outer 
continuum peaks are also real. The different emission distributions are likely due to the combination 
of different amounts of envelope (and, in molecular lines, possibly also ambient) emission recovered, 
differences in optical depths, and molecular abundance variations; these different effects are difficult 
to disentangle with the available data. 

Figures [3] and [4] show velocity channel maps of the 13 CO and C 18 lines, respectively. Around 
Vlsr ~7 km s _1 , both the 13 CO and C 18 emission are significantly suppressed, a likely effect of 
missing flux at the ambient velocity of the molecular gas immediately around L1551 NE. Indeed, 
singl e-dish CS (3-2, 5-4, 7-6) spect ra towards L1551 NE show a central velocity of Vlsr ~7 km 



s _1 ( Moriartv-Schieven et al. 1995bl ). A single-dish C 18 (3-2) spectrum taken towards L1551 NE 
with the CSO shows two distinct components: one having a narrower velocity width (~0.68 km 
s _1 ) with a central velocity of Vlsr ~6.7 km s , and the ot her having a broader velocity width 



(~2.2 km s _1 ) with a central velocity of Vlsr ~7.0 km s _1 ((Fuller fe Ladd 20021 ). The velocity 
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range and the central velocity of the broader component are similar to those of the C 18 emission 
observed with the SMA. We therefore adopt a systemic velocity for L1551 NE of Vlsr = 7 km 
s _1 . In the 13 CO velocity channel maps, there are high-velocity blueshifted (2.5 - 3.6 km s _1 ) 
and redshifted (9.8 - 10.2 km s _1 ) components just displaced to the west and east of Source A, 
respectively. These high-velocity 13 CO components are not detected in the C 18 emission. In 
the lower-blueshifted velocity (~3.9 - 5.8 km s -1 ), both the 13 CO and C 18 emission are located 
predominantly to the north of the protobinary, and in the lower-redshifted velocity (~7.6 km s _1 
- 9.5 km s _1 ) to the south of the protobinary. Furthermore, the extents of the 13 CO and C 18 
emission appear to become progressively larger as the velocity approaches the systemic velocity. 
A comparison of the SMA C 18 (3-2) integrated spectrum with the CSO C 18 (3-2) spectrum 
shows that from ~3.9 to 5.4 km s _1 and from ~8.7 to 9.5 km s _1 almost 100% of the total C 18 
(3-2) flux is recovered in our SMA observations, although the limited signal-to-noise ratio of the 
single-dish spectrum at the line wings prevents us from accurately estimating the amount of the 
recovered flux. Thus, the systematic increase in the spatial extent of the C 18 (3-2) emission over 
these velocity ranges, where the detection level is more than 6a except for that at the bluest and 
the reddest velocities, is probably real. 

To better highlight the spatial-kinematic distributions of the different features traced by the 
13 CO and C 18 emission, we integrated these velocity channel maps over the four different velocity 
ranges discussed above, that is, high-blueshifted (2.5 km s _1 - 3.6 km s -1 ) and redshifted velocities 
(9.8 km s" 1 - 10.2 km s" 1 ) where only 13 CO emission is detected, and low-blueshifted (3.9 km 
s _1 to 5.8 km s _1 ) and redshifted velocities (7.6 km s^ 1 to 9.5 km s" 1 ) where C 18 emission is 
also detected (Figured]). Overall, the low- velocity blueshifted and redshifted components in the 
13 CO and C 18 maps are located to the north-west and south-east, respectively, of the L1551 
NE system. Furthermore, the axis connecting the blueshifted and redshifted emission peaks is 
approximately perpendicular to the axis of the [Fe II] jets driven by Source A (see blue and red 
arrows in Figure [5]). Although the symmetric center of the blueshifted and redshifted emission 
appears to be closer to the position of Source A, pinpointing the precise location of the symmetric 
center may not be straightforward as it seems. A comparison with the CSO C 18 (3-2) spectrum 
indicates that the amount of the missing flux at the blueshifted part of the SMA C 18 emission 
is more than that at the redshifted part (for example, at V sys - 1.0 km s _1 the amount of the 
missing flux is ~80%, while at V sys + 1.0 km s _1 ~60%). Furthermore, the location of L1551 NE 
ov erlaps with the spatially-exten ded redshifted (~7 - ~12 km s _1 ) outflow driven by L1551 IRS 
5 ( Moriarty-Schieven et al. 20061 ). Thus the redshifted part of the 13 CO emission may be more 



severely affected by missing flux than the blueshifted part. These "asymmetric" contributions 
from the extended components may distort the "balance" between the blueshifted and redshifted 
emission distributions. 

In the high- velocity range, the compact (~100 AU) blueshifted and redshifted 13 CO compo- 
nents are located to the north-west and south-east of Source A, respectively. The position angle 
of the axis connecting the high-velocity blueshifted and redshifted 13 CO components is different 
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from that of the low-velocity components; specifically, this axis is not orthogonal to the axis of 
the [Fe II] jets driven by Source A. These results imply that the origin of these high-velocity 13 CO 
components is different from that of the low-velocity components. 

Figure [6] shows intensity-weighted mean- velocity maps of the low-velocity 13 CO and C 18 
emission. As mentioned above, the low-velocity blueshifted emission is distributed north-west and 
redshifted emission south-east of the L1551 NE system. We do not detect a velocity gradient along 
the minor axis of this structure indicative of an inflow, or contamination by an outflow. Thus, both 
the spatial and kinematic major axes of this feature are approximately perpendicular to the axis 
of the associated [Fe II] jets (blue and red arrows in Figure [6j). 



4. Analysis 

As shown in the previous section, the 13 CO (3-2) and C 18 (3-2) emission immediately sur- 
rounding L1551 NE exhibit a conspicuous velocity gradient along their major axes perpendicular to 
the axis of the [Fe II] jets driven by Source A. Such a velocity gradient could be produced by either 
a purely rotating or a rotating and contracting (i.e., infalling component) disk. To differentiate 
between these possibilities, we tried model fittings of geometrically-thin Keplerian or infalling disks 
to the observed 13 CO and the C 18 channel maps. 

Since we assume that the disk is geometrically thin, the centroid line-of-sight velocity at each 
position in right ascension and declination with respect to the disk center (= VLOs( a ,$)) can be 
expressed as; 

VLOs{a, 5) = v sys + v rot (r) cos($ - 9) + v rad (r) sin($ - 9), (2) 

where 

(X 

<J> = arctan(— ), (3) 
o 



:—) 2 + y 2 , (4) 

cost 

x = a cos(#) — 5 sin(0), (5) 

y = asm(9) + £cos(6>). (6) 

In the above expressions v sys is the systemic velocity, 9 is the position angle of the disk major 
axis, and x and y are coordinates of the disk along the minor and major axes, respectively. v ro t(r) 
denotes the rotational velocity of the disk as a function of the radius (= r). In the case of Keplerian 
rotation, v ro t(r) is expressed as; 

v r ot{r) = smiy , (7) 

where i is the inclination angle of the disk from the plane of the sky, G is the gravitational con- 
stant, and M* is the mass of the central star. In the case of rotation with the conserved angular 
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momentum, which is expected for an infalling disk (jOhashi et al. 19971 ; lYen et al. 201ll ). v ro t(r) is 
expressed as; 

j 

Vrot{r) = sini~, (8) 
r 

where j denotes the specific angular momentum of the rotation. In the case of the infalling disk, 
radial motion (= v ra d(r)) is also present and can be expressed as; 



(\ ■ ■ 2GM * fa\ 
Vrad (r) = smtd — - — . (9) 

The velocity channel maps of the model disk (= S mo dei(a,5,v)) can be expressed as; 

Smodei{a,5,v) = (S mom0 (a,5)/aV2TT) x exp(— — — T^pT" 1 ""^ )> ( 10 ) 

where S mom o(a,5) denotes the moment map of the model disk and a the internal velocity dis- 
persion. 

In the model fitting, S rnom o(a, 8) was assumed to be the same as the observed 13 CO or C 18 
moment maps (Figure [2]). We adopted a to be 0.4 km s , inferred approximately from the 
13 CO and C 18 spectra. These assumptions mean that we fit only the global, systematic velocity 
structure in 3-dimensional space, and that we do not fit the 2-dimensional distributions of the 
molecular emission. The center of the model disk was fixed to be the position of Source A, and 
v sys was set to be 7 km s . The 13 CO and C 18 velocity channel maps around the systemic 
velocity (V LS r = 6.52 & 6.89 km s -1 for 13 CO and V LS r = 6.86 & 7.23 km s -1 for C 18 0) were 
excluded from the fitting, since in these velocities the bulk of the emission is "resolved out" . On 
these assumptions, we conducted minimum x 2 -fittings of the model Keplerian disk and the model 
infalling disk to the observed 13 CO and C 18 velocity channel maps, with M*, 9, i, and j (in the 
case of the infalling disk) as fitting parameters; 

x 2 _ y^^ Sobsja, S, v) - S^^ja, 5, v) ^ ^ ^ 

c rms 

a,o,v a,o,v 

where S b s (a,S, v) denotes the observed velocity channel maps, S^*f e f' : ' (a, 5, v) the model velocity 
channel maps with a given M*, 6, i, and j, and a rms is the rms noise level of the observed velocity 
channel maps. 

The middle and lower panels of Figure [7] show the best-fit Keplerian-disk model and the resid- 
ual velocity channel maps of the C 18 emission, respectively, where the best-fit parameters are M+ 
= 0.8^Q4 Mq, 9 = 167 0+ ^ , i = -62 0+ ^ - The residual velocity channel maps do not show any 
systematic features although there are occasional 4a peaks; such random features with 4<7 peaks are 
also seen in the Keplerian model to DM Tau by Guilloteau & Dutrey (1998). The best-fit residual 
rms is ~0.144 Jy beam -1 , which is slightly higher than but comparable to the noise level of the 
observed C 18 velocity channel maps (ler ~0.117 Jy beam -1 ). Our model fit of an infalling disk 
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to the C 18 velocity channel maps, on the other hand, cannot find the local x 2 minimum within 
the range of 6 from 109° to 199°; i.e., the disk major axis is approximately perpendicular to the 
axis of the [Fe II] jets and more than 45° away from the jet axis. These results show that a simple 
geometrically-thin Keplerian-Disk model with only three fitting parameters reproduces the global 
velocity feature in the C 18 velocity channel maps. The position angle of the Keplerian rotation 
is consistent with the major axis of the C 18 total integrated intensity map, and is approximately 
perpendicular to the axis of the [Fe II] jets driven by Source A. The estimated inclination angle 
of the disk from the plane of the sky is also a pproximately perpend icular to the inferred inclina- 



tion angle of the [Fe II] jets of i = 30° - 45° (|Havashi Pvo 20091 ) . The position angle and the 
extent of the Keplerian disk in the C 18 emission are similar to those of the ring-like structure 
in the 335 GHz continuum emission surrounding the binary protostars. These results show that 
the circumbinary continuum feature is indeed a Keplerian disk surrounding the LI 551 NE binary 
system, or "Keplerian circumbinary disk". Interestingly, the estimated position and inclination 
angles of the circumbin ary disk around LI 551 NE ar e similar to those of the circumstellar disk s 



(|Lim Sz Takakuwa 20061 ) and the protostellar envelope (jMomose et al. 19981 ; iTakakuwa et al. 2 004) 
around L1551 IRS 5, but their rotation is in the opposite directions (in L1551 IRS 5 the northern 
part is redshifted). 

By contrast, we could not find a satisfactory Keplerian fit (the peak residual is more than 
10cr) to the 13 CO velocity channel maps. The same Keplerian-disk parameters do reproduce the 
overall spatial-velocity distributions in the 13 CO velocity channel maps, albeit with large residuals, 
presumably because the 13 CO emission is optically thick and likely affected by the contamination 
from the extended envelope component and possibly the outflow from L1551 IRS 5. Our simple 
geometrically-thin disk model, which implicitly assumes optically-thin emission, cannot therefore 
be accurately applied for the 13 CO emission. 

Figure [8] shows Position - Velocity (P-V) diagrams of the 13 CO (left panels) and C 18 emission 
(right) along the major (upper) and minor axes (lower) of the model Keplerian disks, passing 
through the location of Source A. In the P-V diagrams along the major axis, the blueshifted and 
redshifted components are distinctly separated and located to the north and south, respectively, 
of LI 551 NE. The emission at higher velocities is located closer to Source A than that at lower 
velocities, with this trend more pronounced in the 13 CO line. The best-fit Keplerian-rotation 
curve is shown overlaid as solid black curves, and the upper and lower ends of the rotational 
velocity derived from the error bars of the fitting parameters from the x 2 -fitting overlaid as dashed 
black curves. For comparison, a rotation curve with conserved angular momentum as expected 
for infalling gas (j = 6xl0 -4 pc km s _1 ; red curve) is also plotted. Although both black and red 
curves trace the higher-velocity parts close to the center {Vlsr ~2.5 - 4.5 km s _1 & 9 - 10 km 
s -1 ) equally well, the Keplerian-rotation curve better traces the lower- velocity parts (~5 - 6 km 
s _1 & 8 - 9 km s -1 ). We note that the effect of the missing flux may distort the true rotation 
curves. Keplerian-rotation curves (v ro t oc r -0 ' 5 ) are, however, shallower than rotation curves with 
conserved angular momenta (v rot oc r^ 1 ), and thus the spatial extent of the Keplerian disk at a 
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given velocity is larger, as shown by the black and red curves. Hence, it is possible that a genuine 
Keplerian rotation mimics a r _1 rotation due to the effects of missing flux, but it is unlikely that 
the effect of the missing flux distorts a genuine r _1 rotation to a Keplerian rotation. 



5. Discussion 



5.1. Evolution of Disks around Protostars 



Our SMA observations of L1551 NE in the 13 CO and C 18 (3-2) lines have revealed a Keplerian 
disk around the protobinary system. From the outer extent of the molecular-line distributions, 
the outermost detectable radius of the Keplerian disk is estimated to be ~300 AU. The ring-like 
component seen in the 335 GHz continuum emission likely traces this Keplerian circumbinary disk, 
while the central component may arise from circumstellar disks around the binary companions. The 
central mass derived from the Keplerian fitting (~0.8 M ) is likely the mass of the protobinary, 
since the mass of the Keplerian circumbinary disk plus circumstellar disks derived from the 335 
GHz continuum emission (~0.047 M & ) is much smaller. 

Previous interferometric observations of Class protost ars have revealed infal l ing motion in 
compact (5,600 AU) circumstellar structures; i.e., peudo-disks (ITakakuwa et al. 20071 ; iBrinch et al. 2009 : 



Yen et al. 2010 : lYen et al. 20111 ). On the other hand, the large SMA survey for protostellar 



sources ("PROS AC" project) has found Keplerian disks in the HCO + (3-2) line around four 
Class I protostars (L1489-IRS, IRS 63, Elias 29, and IRS 43 ) out of their ten Class I samples 



(jBrinch et al. 20071 : lLommen et al. 20081; I Jprgensen et al. 20091 ). one of which (L1489-IRS) was re 



ported as a protobinary (IHogerheijde &: Sandell 2000 ) . In the other six Class I sources the velocity 
fields traced by the HCO + (3-2) line are rather complicated, presumably due to the contamination 
from the outflows. The reported Keplerian disks have radii ~100 - 200 AU and masses <0.007 to 
0.055 M & , and enclosed stellar masses 0.37 - 2.5 M & . These results imply that Keplerian disks 
can be present around Class I protostars as well as T-Tauri stars (e.g., Guilloteau & Dutrey 1998; 
Guilloteau et al. 1999). The measured physical properties of these protostars and the Keplerian 
disks, including those of L1551 NE, are summarize in Table [2j 

Table [2] shows that there are no clear differences of the disk radius and the enclosed stellar 
mass between L1551 NE and the other sources. On the other hand, the bolometric temperature 
of L1551 NE (Tboi = 91 K) is a factor >3 lower than that of the other protostars, and the mass 
of the surrounding envelope around L1551 NE is a factor of ~4 higher than that of L1489-IRS 
and over an order of magnitude higher than that around the other protostars. Furthermore, the 
fraction of the central protostellar mass to the total (star + disk + envelope) mass is the lowest 
in L1551 NE (= 65%). These results imply that either L1551 NE will evolve into a much more 
massive star compared with the other protostars listed in Table [2j or L1551 NE is the youngest 
Class I pro tostar yet found to have a Keplerian disk . Since the L1551 region forms solar-type (< 1 
M ) stars (jHavashi et al. 19931 : llVIomose et al. 19981 : ITakakuwa et al. 20041 : lHavashi fc Pvo 20091 ). 
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the former case is unlikely. Thus, we suggest that LI 551 NE has experienced "transition" from 
an infalling, pseudo-disk to a Keplerian disk most recently among the Class I sources with the 
Keplerian disks. The formation mechanism of Keplerian disks around protostars will be discussed 
in the next section. 



5.2. On the Formation of a Keplerian Circumbinary Disk 

Our SMA observations of L1551 NE provide firm evidence that Keplerian circumbinary disks 
are present in addition to Keplerian circumstellar disks around Class I protostars. On the theoretical 
side, however, it is still difficult to form Keplerian disks around either single or multiple protostars. 
Mellon & Li (2008; 2009) have conducted 2-D axisymmetric simulations of the collapse of rotating, 
magnetized singular isothermal cores. They found that, even in the case of weak initial magnetic 
fields (mass-to-flux ratio A <100), magnetic fields efficiently carry angular momenta outwards, 
preventing the formation of a rotationally-supported disk. They showed that a field strength in the 
dense core L1544 inferred observationally (A <10; Crutcher & Troland 2000) completely suppresses 
the formation of the centrifugally-supported disk, even if the ambipolar diffusion is taken into 
account. They concluded that additional processes, such as Ohmic dissipation, Hall effect, or 
the dispersion of the surrounding envelope that anchors the magnetic brake, are required for a 
Keplerian disk to form. Machida et al. (2010; 2011a; b) have performed 3-D magnetohydro dynamic 
simulations of a Bonnor-Ebert sphere with a mass of 1 Mq and different initial rotations and 
magnetic fields, and have calculated the evolutions. Their simulations show that the first adiabatic 
core formed prior to the formation of the protostar becomes a centrifugally-supported disk with 
a size of >100 AU by the end of the calculations. Their results show that, however, the mass of 
the formed disk is considerably (~factor 2 - 100) larger than that of the central protostar (see 
also Vorobyov 2009), and hence the formed disk is subject to further fragmentation. Further long- 
term simulations, including calculations of fragmentations of the massive disks and subsequent 
evolutions, are required to theoretically reproduce the observed Keplerian disks around the Class I 
protostars. 

Although a Keplerian circumbinary disk is identified ar ound L1551 NE, L1551 NE is associated 



with active jets (jReipurth et al. 2000i : lHavashi fc Pvo 20091 ) and molecular outflows (|Moriartv-Schieven et al. 1995 



suggesting presence of accreting circumstellar disks around the binary companions. The central 
compact dust emission (Figure [T]) may indeed arise from the circumstellar disks. Simulations of 
fragmentations of a first core or a massive disk into multiple protostellar systems (i.e., Machida 
et al. 2008) and the further evolutions are required to reproduce a Keplerian circumbinary disk 
plus accreting circumstellar disks in a protobinary system. In the case of binary protostars, tidal 
effects from the two protostars can transport rotational angular mome nta outwards and create a 



'gap" in the circumbinary disk (IBate Bonnell 19971 ; lOchi et al. 20051 ). Furthermore, the inward 



transportation of magnetic fields in binary systems may be quite different from that around single 
protostars, because the mass accretion toward the binary protostars can occur through "spiral" 
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accretion str eams that connects between the inner e dge of the circumbinary disk and the circum- 
stellar disks ( Bate fc Bonnell 19971 ; lOchi et al. 20051 ) . The high-velocity 13 CO components shown 



in Figure [5] may correspond to such accretion streams. If the angular momentum transported out- 
wards piles up at the inner edge of the circumbinary disk, this disk can then assume Keplerian 
rotation. 

Finally, we will comment on an interesting application of Keplerian circumbinary disks. By 
measuring the position of the Keplerian rotation center with respect to the binary positions, in 
principle it is also possible to derive binary mass ratios, orbital radii, and orbital periods of binary 
protostars. Our x 2 model fitting of the geometrically-thin Keplerian disk with the fixed rotational 
center at the position of Source A provides a reasonable fitting result. This result may suggest 
that the mass of Source A is higher than that of Source B, although our search for the minimum 
X 2 point of the rotational center failed to obtain a statistically significant result (at the positions 
of Sources A and B the reduced x 2 values are ~1.5 and ~1.8, respectively). If we could pinpoint 
the rotational center, the binary mass ratio can simply be estimated from the ratio of the distances 
from the rotational center to the primary and to the secondary. This is potentially a strong method 
to measure mass ratios of protostellar binaries, since previous multi-epoch (> 5 yr), high-resolution 
(^Cf.'l) observations of proper motions of protobinary systems could only deduce the reduced masses 



(jLoinard et al. 20021 ; iRodriguez et al. 20031 ; iLim k, Takakuwa 20061 ) . We can also measure the de- 



projected orbital radii assuming that the binary stars are in the same plane of the circumbinary 
disks. From the binary masses, mass ratios, and the de-projected orbital radii the orbital periods 
can also be obtained. These measurements could put stringent constraints on theories of binary 
formation. 



6. Summary 

We have carried out sub-arcsecond resolution SMA observations of the protobinary system 
L1551 NE in the 335 GHz continuum, 13 CO (3-2), and the C 18 (3-2) emission. The main results 
are summarized below: 

1. The 335 GHz continuum emission in L1551 NE shows an intense central peak close to 
the protobinary, and secondary peaks to the north and south of the protobinary system located 
almost symmetrically with respect to the position of the protobinary. The central component 
may arise from circumstellar disks around each binary companion. The northern and southern 
components comprise a ~600 AU scale circumbinary disk, approximately perpendicular to the axis 
of the associated [Fe II] jets driven by Source A. The mass of the circumstellar disk is estimated to 
be ~0.02 - 0.09 M , and that of the circumbinary disk ~0.03 - 0.12 M . 

2. The extent and position angle of the 13 CO (3-2) and C 18 (3-2) emission are similar to 
those of the 335 GHz continuum emission, although the locations of the emission peaks are different 
among the 13 CO, C 18 0, and the 335 GHz emission. Both the 13 CO and C 18 emission show a 
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systematic velocity gradient along the major axis, which is perpendicular to the axis of the [Fe II] 
jets, such that their blueshifted and the redshifted emission are located to the north and south of the 
protostellar binary, respectively. In addition to this global velocity gradient, the high-velocity (>3 
km s _1 ) 13 CO emission exhibits compact (~100 AU) blueshifted and redshifted components to the 
north-west and south-east of Source A with a different position angle from that of the circumbinary 
disk. We speculate that the high-velocity 13 CO components correspond to accretion streams that 
connects between the inner edge of the circumbinary disk and the circumstellar disks. 

3. x 2 model fitting of geometrically-thin Keplerian disks to the observed 13 CO and C 18 
velocity channel maps was conducted. The C 18 velocity channel maps are well-reproduced by 
a geometrically-thin Keplerian disk model with a central stellar mass of 0.8±g;| M Q , disk posi- 
tion angle of 167 27 o, and a disk inclination angle of -62 17 o, and a rotational center at the 
position of Source A. On the other hand, the 13 CO velocity channel maps cannot be fitted with 
any geometrically-thin Keplerian disk model. This is probably because the 13 CO emission is opti- 
cally thick and is likely affected by the contamination from the outflow and the extended envelope 
component. 

4. Among the Class I sources with identified Keplerian disks, L1551 NE is likely the youngest 
source, since L1551 NE has the lowest bolometric temperature (~91 K), highest envelope mass 
(~0.39 M & ), and the lowest ratio in stellar mass to envelope + disk + stellar mass (~0.65). 
We suggest that L1551 NE has just passed from the "pseudo-disk" phase seen around Class 
protostars to the Keplerian-disk phase. Our SMA observations of L1551 NE provide firm evidence 
that Keplerian disks around binary protostars, Keplerian circumbinary disks, are present as well 
as Keplerian disks around single protostars. 

5. Current theoretical studies cannot reproduce observationally-identified Keplerian disks 
around single or multiple protostars. In the case of binary protostars, tidal effects from two proto- 
stars can transport rotational angular momenta outward, and pile up the angular momenta at the 
inner edge of the circumbinary disk effectively. Furthermore, the inward transportation of mag- 
netic fields in binary systems may be quite different from that around single protostars, because 
the mass accretion toward the binary protostars can occur through "spiral" accretion streams from 
the inner edge of the circumbinary disk to the circumstellar disks. These mechanisms could create 
a Keplerian circumbinary disk around binary protostars. Higher-resolution observations of L1551 
NE with ALMA and / or EVLA, which can separate between the circumstellar disks and the cir- 
cumbinary disk and resolve the accretion streams connecting between the circumbinary disk and 
the circumstellar disks, are required to fully understand the origin of the Keplerian circumbinary 
disk plus the circumstellar disks. 

We are grateful to P. T. P. Ho and N. Ohashi for their fruitful discussions. We would like to 
thank all the SMA staff supporting this work. S.T. acknowledges a grant from the National Science 
Council of Taiwan (NSC 99-2112-M-001-009-MY3) in support of this work. 
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Right ascension (J2000.0) 

Fig. 1. — 335-GHz continuum image of L1551 NE observed with the SMA. Contour levels are from 
3a in steps of 4a until 35cr, and then in steps of 20a (la = 3 mJy beam -1 ). The highest contour 
level is 115<t. Crosses indicate the positions of Source A and Source B as labeled. A filled ellipse 
at the bottom-right corner shows the synthesized beam ((X'80 x 0".b4; P. A. = -87°). Blue and red 
dashed curves delineate the tilted [/-shaped features, and blue and red arrows show the direction 
of the blueshifted and redshifted [Fe II] jets driven by Source A. 
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Fig. 2.— Moment maps of the 13 CO (3-2) (Left) and C 18 (3-2) lines (Right) in L1551 NE 
observed with the SMA. In the 13 CO map, contour levels are from 2a in steps of 3a (la = 0.189 
Jy beam -1 km s™ 1 ), and the integrated-velocity range is 2.28 - 10.76 km s" 1 . In the C 18 map, 
contour levels are from 2a in steps of 2a (la = 0.173 Jy beam -1 km s _1 ), and the integrated- 
velocity range is 3.72 - 9.63 km s _1 . Crosses indicate the positions of the protobinary, and filled 
ellipses at the bottom-right corner the synthesized beam (O^S x / /66; P.A. = -88°). Blue and 
red arrows show the direction of the blueshifted and redshifted [Fe II] jets driven by Source A. A 
righthand scale-bar shows the intensities of the gray-scale images in Jy beam -1 km s . 



-20- 



1 



I 1 1 1 I 1 (' f ' 



V 



1 



3p 



i k , rrr: 



03 



»■ c> 




C0_ 
in ° 



3 - 





O 

CO 



10 





J *0 



A 





3 





% 



o 



0) 



CO 





l\,,ICJ 



- N 
CD 



"J 

/) 



cr in 



03 
O 



o 



- r- 

CD 



1 1 

" : ; 
J . ,1 




a> 



- CO 
CO 



4l 



P 3 



_ 

o 



o 



*0 



05 
O 



•II 







CO 
- iCO 
CO 







o 



(OOOSf) 33.Q 



(ooo2f) oaa 



(ooosr) Daa 



(ooosr) oaa 



Fig. 3 — Velocity channel maps of the 13 CO (3-2) line in L1551 NE observed with the SMA. 
Contour levels are from 2a in steps of 2a (la = 0.107 Jy beam -1 ). Crosses indicate the positions 
of the protobinary, and a filled ellipse at the bottom-right corner in each panel the synthesized 
beam (0 / /95 x / .'66; P. A. = -88°). A number at the top-left corner in each panel denotes the LSR 
velocity. 
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Fig. 4. — Velocity channel maps of the C 18 (3-2) line in L1551 NE observed with the SMA. 
Contour levels are from 2a in steps of 2a (la = 0.117 Jy beam -1 ). Crosses indicate the positions 
of the protobinary, and a filled ellipse at the bottom-right corner in each panel the synthesized 
beam ((Z/95 x C//66; P. A. = -88°). A number at the top-left corner in each panel denotes the LSR 
velocity. 
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Fig. 5. — a) Distributions of the low-velocity blueshifted (blue contours; 3.72 - 5.93 km s -1 ) and 
redshifted C 18 (3-2) emission (red contours; 7.41 - 9.63 km s -1 ) in L1551 NE observed with 
the SMA. Contour levels are from 0.212 Jy beam" 1 km s -1 in steps of 0.424 Jy beam -1 km s -1 
(la = 0.106 Jy beam" 1 km s -1 ). Crosses indicate the positions of the protobinary, and a filled 
ellipse at the bottom-right corner the synthesized beam (0l!95 x / .'66; P.A. = -88°). Blue and 
red arrows show the direction of the blueshifted and redshifted [Fe II] jets driven by Source A. b) 
Distributions of the low- velocity blueshifted (blue contours) and redshifted 13 CO (3-2) emission 
(red contours) in L1551 NE observed with the SMA. The integrated velocity ranges, contour levels, 
and the symbols are the same as those in a), while the la rms noise level is 0.097 Jy beam -1 
km s . c) Distributions of the high- velocity blueshifted (blue contours; 2.24 - 3.72 km s -1 ) and 
redshifted 13 CO (3-2) emission (red contours; 9.63 - 10.37 km s -1 ) in L1551 NE observed with the 
SMA. Contour levels are from 0.158 Jy beam -1 km s -1 in steps of 0.158 Jy beam -1 km s -1 , where 
the la rms noise levels are 0.079 Jy beam -1 km s -1 and 0.056 Jy beam -1 km s -1 in the blueshifted 
and redshifted velocity ranges, respectively. 
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Fig. 6.— Moment 1 maps of the 13 CO (3-2) (Left) and C 18 (3-2) lines (Right) in L1551 NE 
observed with the SMA. Contour levels are in steps of 0.4 km s _1 . The bluest and reddest contour 
levels to the north-west and south-east of Source A in the 13 CO map are 5.2 km s _1 and 8.8 km s _1 , 
respectively, and those in the C 18 map 5.2 km s _1 and 8.4 km s _1 . Crosses indicate the positions 
of the protobinary, and filled ellipses at the bottom-right corners the synthesized beam (Cf.'95 x 
C'66; P. A. = -88°). Blue and red arrows show the direction of the blueshifted and redshifted [Fe 
II] jets driven by Source A. 
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4 h 31 m 44 s 7 44 s 3 
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Fig. 7. — Result of the x 2 fitting of the geometrically-thin Keplerian Disk to the C 18 (3-2) 
velocity channel maps in L1551 NE at the blueshifted (blue contours) and redshifted (red contours) 
velocities. Upper, middle, and lower panels show the observed, model, and the residual velocity 
channel maps, respectively, where the best fit parameters are M* = 0.8 M , 6 = 167°, and % = 
-62°. Contour levels are from 2cr in steps of 2cr (la = 0.117 Jy beam -1 ). Crosses show the position 
of the protobinary, and the filled ellipses at the bottom-right corners the SMA synthesized beam 
(C/95 x 0'/66; P. A. = -88°). 
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Fig. 8.— Position - Velocity (P-V) diagrams of the 13 CO (3-2) (Left) and C 18 (3-2) lines (Right) 
along the major (P.A. = -13°; Upper) and minor axes (P.A. = 77°; Lower) of the best-fit model 
Keplerian disk, passing through the position of Source A. Contour levels are from 2a in steps of 2a 
(la = 0.107 Jy beam -1 and 0.117 Jy beam -1 for the 13 CO and C 18 P-V diagrams, respectively). 
Horizontal and vertical dashed lines denote the position of Source A and systemic velocity of 7 km 
s , respectively. Solid and dashed black curves show the best-fit Keplerian rotation curve and the 
upper and lower ends of the rotation curve derived from the error bars of the fitting parameters, 
while red curves denote rotation with the conserved specific angular momentum expected in infalling 
gas (see text). 
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Table 1. Parameters for the SMA Observations of L1551 NE 



Parameter 


Value 




2010 January 15 2010 January 18 


Number of Antennas 


8 7 


Right ascension (J2000.0) 


g 4 h 31 m 44 s 47 


Declination (J2000.0) 


18° 08' 32'.' 2 


Primary Beam HPBW 


~37" 


Synthesized Beam HPBW (Line) 


0'.'95 x Of.'ffi (P. A. = -88°) 


Synthesized Beam HPBW (Continuum) 


Of.' 80 x 0^'54 (PA. = -87°) 


Baseline Coverage 


29 - 249 (kA) 


Conversion Factor (C 18 0) 


1 (Jy beam- 1 ) = 17.9 (K) 


Frequency Resolution 


406.25 kHz ~0.37 km s" 1 


Bandwidth 


7.74 GHz 


Flux Calibrator 


Titan Callisto 


Gain Calibrator 


0423-013, 3cl20 


Flux (0423-013 Upper) 


3.46 Jy 3.62 Jy 


Flux (0423-013 Lower) 


3.53 Jy 4.42 Jy 


Flux (3cl20 Upper) 


0.80 Jy 0.92 Jy 


Flux (3cl20 Lower) 


0.83 Jy 1.08 Jy 


Passband Calibrator 


3c273 


System Temperature (DSB) 


~300 - 700 K ~400 - 1000 K 


rms noise level (Continuum) 


3 mjy beam 


rms noise level (C ls O) 


0.117 Jy beam" 1 


rms noise level ( 13 CO) 


0.107 Jy beam" 1 



Table 2. Protostellar Sources with Keplerian Disks 



Source 


Lbol 
(L e ) 


Tbol 
(K) 


(AU) 


(Mq) 


(Mq) 


M d 
(Mq) 


M star 

M st ar+M disk + M C nv 
(%) 


references 


L1551 NE 


4.2 


91 


300 


0.8 


0.026 f 


0.39 


65 


1,2,3 


L1489-IRS 


3.7 


238 


200 


1.35 


0.004 


0.093 


93 


4,5 


IRS 43 


6.0 


310 


140 


1.0 


0.0081 


0.026 


97 


5 


IRS 63 


0.79 


351 


100 


0.37 


0.055 


0.022 


83 


5,6 


Elias 29 


13.6 


391 


200 


2.5 


<0.007 


0.025 


>99 


5,6 



a Outer radius of the Keplerian disk. 

b Central protostellar mass derived from the Keplcrian-Disk model fitting. 
c Mass of the Keplerian disk measured from the SMA continuum flux. 

d Mass of the protostellar envelope estimated from the Single-dish sub-millimeter continuum flux. 

References: (1) This work; (2) Frocbrich 2005; (3) Moriarty-Schicvcn et al. 2006; (4) Brinch et al. 2007; (5) 
j0rgensen et al. 2009; (6) Lommen et al. 2008 

f Excluding the mass of the central compact component that may arise from the circumstcllar disks around the 
binary companions. 



